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ABSTRACT 

A model of a TWS radar is developed that provides a realistic 
computer simulation for comparing various radar tracking methods. 

Prediction accuracy of a simplified alpha - beta tracker is 
compared to that of an adaptive filter. In addition, the effect 
on radar tracking of a variable gate size correlation technique 
is investigated. 
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CHAPTER I 



INTRODUCTION 

A need exists for more extensive use of computer simulations 
as an aid for defining and obtaining the solutions to systems engin- 
eering problems. To gain some insight and experience in this field 
a computer simulation of a hypothetical point defense radar system 
is carried out. 

The criteria for defining the radar parameters and method of 
operation were derived from the following situation. 

A destroyer class ship should be able to' provide itself with 
adequate protection against a multiple aircraft raid while transiting 
singly. In order to provide timely detection and adequate intelligence 
for launch of a defensive weapon the shipboard radar must be capable 
of the following functions: 

1. Full hemispherical coverage. 

2. Medium range search i.e. 50-100 n.m. 

3. Multiple target track while scan. 

High .data rates in critical threat sectors. 

In addition, to minimize costs the radar design should be compa- 
tible with existing weapons. To provide a realistic approach, only 
those designs and techniques capable of being implemented with curren- 
tly available components were allowed. 

A phased array antenna with a monopulse beam was chosen to meet 
the requirement of high data rates on multiple targets. Phased arrays 
are inherently limited in scan. To meet the requirement of full 
360 ® coverage it was necessary to rotate the antenna mechanically in 
azimuth. Formulation of the problem in this way led to the selection 
of two modes of operation. The first is characterized by complete 
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radar coverage with mechanical rotation in azimuth and electronic 
scanning in elevation, with monopulse in elevation angle only. This 
mode has limited tracking capability due to a relatively low revisit 
rate. Mode 2, a sector scan with an electronically steered beam in 
azimuth and elevation and monopulsing in two coordinates, is capable 
of providing 0° to 90® elevation coverage and ^ ^5® in relative az- 
imuth . 

It is felt that the two-mode solution is justified since one can 
assume with some confidence that the raid will be radially inbound 
with little or no intentional maneuvering at ranges greater than twenty- 
five nautical miles. Once the raid has been detected and evaluated as 
such by the observer. Mode 2 can be selected to provide a high data 
rate to track the target, whose probability of maneuver increases as 
the range to the target decreases and as his weapon launch point is 
approached. 

The switching range is not critical and can be assumed to be a 
function of the maximum range of the expected weapon available to 
the enemy aircraft (of the order of 5 to 15 n.m. ) plus some safety 
time. 

^ Once the radar was modeled and the computer program written, the 
simulation became a convenient vehicle for investigation of the radar 
tracking equations which are a critical facet of a T W S system. A 
comparison of prediction accuracy and resolution of target tracks 
between a classical alpha/beta tracker and an adaptive technique was 
pursued as an extension of the basic premise of utilizing the computer 
as a device for design and analysis. 
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CHAPTER II 



RADAR MODEL 



A discussion of some assumptions which affected the choice of 
the type of radar to be modeled was presented in the introduction. 

This chapter presents some of the reasoning and methods used to select 
basic parameters for the radar necessary to establish a reasonable 
time scale for the simulation. This initial section is followed 
by a discussion of a few of the more important subroutines in enough 
detail to supplement the flow diagrams of appendix II where necessary. 
Choice of Basic Parameters. 

Assume: (l) Target velocity = 600n.m./hr. 



( 2 ) Attack profile shown in fig. 

(3) Target weapon release point lO m.m. max. 

(h) Defensive weapon - missile (mach 3). 

(5) Missile requires a high data rate for command guidance. 

( 6 ) Antenna height (h^^) = 50 ft. 



ALT (feet) 



.. 2500 



2000 



TARGET ASSUMED TO BE ON 
A RADIALLY INBD. COURSE. 








5 



10 



15 

(Nautical Miles) 



20 



25 



Fig. 1 ATTACK PROFILE 
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Desiring a 30 second sector scan in which to fire and control 



the defensive missle, we must commence sector scan at a minimum range 
of 20 n. m. This ensures an impact point of approximately 15 n . m. , 
providing adequate time for a second attempt before the attacking 
aircraft reaches his weapon launch point. Arbitrarily assuming twice 
this range (40 n. m. ) will provide adequate time for threat evaluation. 
Based on the radar horizon given by the following equation, 50 n . m. 
was chosen as the range at which a probability of detection, equal 
to 90^ on a single scan, would provide sufficient initial detection 
capability. 



Having obtained an idealized Rmax = n. m. , we can calculate 
the pulse repetition frequency. 



(miles ) 

h^ = height of target A/C. 
given : h^ = 50 ft. 




r Ul miles h^ = 500 ft. 
r ^ 5^ miles h^ = 1000 ft. 



R(unamb.) - c = ^5 n. m. 



2fr 

f^ = 3 x1qQ m/sec ^ 5-^ X 10 
2x h3 n.m. 



n.m. 



m. 



fj, = l800 cycles/sec. 



12 



Assiaming a fan shaped antenna beam (l^ azimuth by 5.5° in 
elevation) and requiring the equivalent of integrating h pulses per 
scan we determine the time for one revolution of the antenna in 
mode 1 . 



1 _ h cycles 

1800 cycles/sec. 1 ° Beam 

Position 



360^ 

REV. 



0.8 SEC . 
REV. 



This figure must be multiplied by the number of 5.5° sections in 
elevation that must be covered. Taking into account the target 
profile and realizing the advantages of short revisit rates a 
compromise value of 5^.5° coverage in elevation requiring 8 sec. 
total scan time was' selected. 
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To simulate the continuous motion of the antenna, azimuth is 
divided into 0.02° increments. The following diagram describes 
the antenna motion. 




Fig. 2. SIMULATED ANTENNA MOTION — MODE 1 



Turning now to a discussion of mode 2 or full electronic sector 



scan the following assumptions were made: 

(1) Defensive missile has a 20 millisecond command guidance 
requirement repeated at intervals of 1 second. 

(2) Special purpose heam steering computer capable of: 

(a) Calculating 200 beam positions per second. 

(b) superimposing I collimation commands on the beam 
position to obtain variable beam shape. 

(c) Beam shape = 2° x 2® pencil beam. 

(3) 5 ms. illumination time per beam position except 25 ms. 
illumination time for beam positions where targets are predicted 
to be present. 

(4) 10 - target capability > maximum. 

The above assumptions lead to the establishment of the following 
table: 

TABLE 1. Search Beam Positions vs. known targets 
Targets 

0 
2 

6 
8 

10 



Search Beam Positions Remaining Per 
1 Second Interval. 

200 

190 

l80 

170 

i6o 

150 
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In a i 45® AZ. by 60® EL, sector we have 45 x 30 = 1350 (2® x 2®) beam 

positions. For the worst case of 10 targets engaged the revisit 

rate to a "same” search position would be 13^0 = approximately 9 

150 

seconds. This scan interval is reduced to 6.75 seconds for the case 
of 0 targets engaged. 




Fig. 3. BEAM SCANNING SEQUENCE FOR 2 TARGETS ENGAGED. MODE 2. 



The target detection model for mode I and mode II are essentially 
identical and can be described in the following way. If the radar 
beam passes over a target the target range as given by the target 
generator is compared to the values in Table 2 to determine the 
probability of detection for that range. If Random, the result of 
calling the uniform random number generator RANI, is less than the 
probability of detection, we assume a detection; otherwise a miss 
or noise return. 
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TABLE 2. Probability of detection vs. range. 



Prob. of Det. Target Range 



*: * 



0.22 


6o.O 


0.5 


55.0 


0 . 6 


5it.O 


O.T 


53.0 


0.8 


52.0 


0.85 


51.0 


0.9 


50.0 


0.95 


49.0 


0.99 


48.0 


0.999 


4?.0 



Having established a radar target the radar range RADRNG is 

‘•I 

computed by the equation RADRNG = TARGRNG + RANERR where RANERR is 
obtained by calling the GAUSSIAN random number generator and using 
the equation: 

RANERR = DEV * 0.3 

We can see from the above equations that the radar range is 
equal to the true range plus an error term with a standard deviation 
by choice, of 0.3 n. m. 

It should be noted at this point that all radar measurements 
are computed in essentially the same manner. That is RADRNG, AZRAD 
and ELRAD are computed as noisy observations of the ”true” target 
position, as given by the target generation subroutine, TARGEN. 
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CORRELATION: 



To correlate current radar observations with tracks established 
from predictions on previous observations, subroutine corrass was 
devised. Current radar observations are stored in the matrices RADRNG, 
AZRAD, ELRAD^ While predicted values from the tracker are stored 
similiarly. Once during each revolution of the antenna each track 
is compared with each observation; for each observation that a track 
agrees with, a one is placed in the corresponding location of the 
C(l,j) correlation matrix. The columns of the correlation matrix 
correspond to tracks and the rows correspond to observations. For 
a track to agree with an observation certain tolerances were assumed: 
range must be within t X n.m., azimuth and elevation angles must be 
within t Y degrees . 



X and Y are determined as follows: 
(l) Alpha - Beta filter. 



Range gate 
Azimuth gate 
Elevation gate 
(2) Kalman Filter. 
Range gate 
Azimuth gate 
Elevation gate 



= 2.5 n.m. 

= 2.5 n.m. 

= 1.5 n.m. 

= 2.0 * SQRTF lP(l,l)] n.m. 
= 1.0 * SQRTF [P(3,3)] n.m. 
= 2.5 * SQRTF [P(5,5)] n.m. 



After this initial correlation several observations may correlate 
with one track or several tracks may correlate with one observation. 

In these situations, a set of four association rules are applied to 
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the correlation matrix, C (l,j), to assign observations to the 
proper tracks . 

These rules implemented by subroutine ASSOC, are as follows: 

(1) A track which correlates with several observations rejects 
any observation held in common with another track, if the common 
observation is the only observation correlating with the second track. 

(2) A track correlating with several observations, some of 

V 

which are not held in common with other tracks, rejects observations 
which are held in common with other tracks. 

(3) When several observations correlate with one track, the 
closest observation is associated with that track. 

(U) When several tracks correlate with one observation, the 
observation is associated with the closest track. 

Rules 3 and h are based on the value of range since it has the 
smallest for measurement noise. The above correlation and assoc- 
iation routine is essentially the same as that given by [T]. 
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CHAPTER III 



FILTERS 



Smoothing of raw radar reports and prediction of future obser- 
vations, basic functions of an automatic tracking system, can be 
performed conveniently by a set of equations implemented on a com- 
puter. In this simulation two distinct schemes were chosen to 
perform these functions. In the simplest case target tracks are 
based on smoothing and prediction of an alpha-beta tracker operating 
in a cartesian coordinate reference frame. 

The smoothing equation is: 



/\sn ~/^pn ■'■(X 0\n 



The prediction equations : 



An ~An-l ■^•^^An "Apn^ 

m 

y\_pn+l "Asn An * 



Where : 




Smooth value of X,Y or Z for the n*^^ scarf. 




. Predicted value of X,Y or Z for the scan. 

^ = Noisy observation of X,Y or Z for the n"^^ scan. 

= Smooth prediction of X,Y or Z component of velocity 



scan. 



for the n*^^ 



scan. 
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Where : 



A = Predicted X,Y or Z coordinate for the N + 1 scan. 

A Pn+l 

T = Time between looks which is essentially constant. 

T =1 second for sector scan. 

T =8 seconds for full scan. 

= The smoothing parameter. 

P 

(3 = <x for optimum filtering. 

2 - 0 ( 

Normally associated with the alpha-beta smoothing prediction 
equations are a set of rules for determining bin size. The rules 
are usually based on track firmness, length of time since last 
correlation and radar range. Basically this set of rules for 
varying the correlation gate size is an attempt to make the system 
adaptive to target dynamics. 

Experimentation with the simple alpha-beta tracker indicated 
that given a reasonable fixed gate size based only on accuracy of the 
radar measurements, tracks could be maintained on targets with mild 
maneuvers, ie., less than 3^/sec. turns, with little difficulty given 
a suitable value of alpha. An important result observed from these 
simulations was that the real problem associated with tracking is 
the time delay before sensing a target’s maneuver initiation. 

Investigation of a criterion for determining target maneuvers 
seemed to be the next logical step. Changes in the magnitude and 
angle of the target velocity vector in the X - Y plane were looked 
at briefly as one possibility. From the data accumulated, inciden- 
tal to other tasks being performed, the changes in angle appeared 
to be an unlikely candidate due to extreme noise caused by measure- 
ment errors of the radar. The magnitude of the vector, on the other 
hand, appeared to be rather insensitive to target maneuvers. 
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In lieu of pursuing this somewhat uncertain investigation, it was 
decided to use techniques that were in comparison well defined. The 
Kalman filter was chosen as the best solution to the problem for the 
following reasons: 

(1) While the Kalman filter can not predict target maneuvers, 
the P matrix generated by the filter equations is directly influenced 
by target maneuverability (Q matrix) and at the same time is a 
measure of our confidence in predicted values. It seemed reasonable, 
therefore, to make the correlation gate size proportional to the 
corresponding values of the P matrix. 

(2) If an observation is missed due to a target maneuver we 
can expand the gate size to some larger value to increase the pro- 
bability of correlation on the next scan with automatic gate size 
reduction if we succeed. 

The Kalman Filter Equations as given in (U): 



" Pn/n-1 

= ^/n-1 * <=n 

Bn/n ” Bn/n-1 - '^n^Bn/n-l 
t*l/n 

Bntl/n = * « 
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These equations represented in block diagram form appear in 



fig. 




Fig. BLOCK DIAGRAM OF FILTER EQUATIONS. 






where : 



or 



Z = H X + V 



Z = 



Radar Range 
Radar Azimuth 
Radar Elevation. 



V— ^ represents the noise associated with the radar measurement 



of the three observables and is gaussian. 



or 



H = observability matrix. 



H = 



1 0 0 0 0 0 
0 0 1 0 0 0 
0 0 0 0 1 0 



0 = The state transition matrix that represents the dynamics 



of the aircraft target. Assuming the aircraft can be rep- 



1 



resented as pure inertia of a -q 



plant . 
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We have 



u(sj ^ 



1/s 



X2(0) 



1/S 



1/S. 



Xj^Co) 

1/s 



X5(S) 



X2(S) 



Xt(S) 



X]^(S) 



Yields : 



Fig. 5 . FLOW GRAPH OF l/S PLANT. 



X^(S) = 0 X^ (S) + 1 X 2 (S) + 0 U(S) 



X^(S) = 0 X^ (S) + 0 X 2 (S) + 1 U(S) 



hence; 



0 1 - 

k ^ ^ [ si-a ] = 

0 0. 



S -1 

0 1/s 



where : ^ A = 

-1 P °1 " F 

[-A] =Li sJ = Lo sJ = 



[si- 



A 



i/s i/s^ 

0 l/S 



we know: ^ (S) = [SI-A] ^ 

^ t1 



1/S 1/S^ 

0 1/s 



hence: 0 (t) = 



0 1 



C(S) 



Assuming the components of the X vector are uncoupled. 



0(t) 



1 


t 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 


t 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 


t 


0 


0 


0 


0 


0 


1 



R = The covariance of the measurement aoise is by defini- 
tion = E [V V^] 



In this case the- standard deviation of the measurement 
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noise on Range, Azimuth and elevation are known to have 



the following values. 
Range = 0.3n.m. 

Az =0.5 deg. 

El =0.5 deg. 

hence : 

R = 



0.09 0 0 
0 0.25 0 
0 0 0.25 



n/n~l represents the predicted states of a radar target. 

To initialize the filter for a particular target, that is specify 

X 1/0, the following method is used. 

Radar Range 
R = -0.166 



X 1/0 



Radar Azimuth 
AZ = 0.0 



Radar Elevation 
El = 0.0 



This method assumes the target is closing the radar radially 
(az = O.O) at a low altitude (El’^ O.O) with a velocity of approx- 
imately 600 n.m. /hr. 

P n/n-1 represents the amount of uncertainty in the predicted 
values of target states X n/n - 1. 

P = E [(X - , ) (X - X / , )'^] 

where: X is the true state. 



X is the predicted value. 

To initialize the filter we must provide. 

^1/6 “1/0 previously and R the covariance of the meas- 

urement noise; 



P (l,l), P ( 3 , 3)5 P ( 595 ) are known. 

Again relying on the assumption of a radially closing target at low 
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altitude with a speed of 600 rx.m./hr. we can assume the remain- 

2 2 2 

ing diagonal elements of P which represent 

will he small. Actual values were determined empirically and as a 



result the initial covariance of error matrix is determined. 



P = 



0.09 0 0 0 0 0 

0 0.0278 0 000 
0 0 2.5 0 0 0 

0 0 0 0.25 0 0 
0 0 0 0 0.25 0 
0 0000 0.00109 

It should be noted here that = 2.5 is used in place 

of 0.25 to aid tracking maneuvering targets. If an observation fails 
to correlate with an established track the P matrix for that track 
is set equal to the above values. The gate sizes, and in particular, 
the Azimuth gate, are directly related to the P matrix values. This 



procedure enlarges the gate size thus increasing the probability of 
correlation on the next scan. 

Lastly but certainly not of the least importance is the Q matrix 
specification. The Q matrix which represents perturbations to the 
state vector due to target maneuvers must be representative of the 
full range of movements the target is expected to use. 

Values for Q were estimated by assuming’ a ”worst” case such as 
6000 ft. /min. rate of climb, 1 g. linear acceleration and 12° / 
second turning rates. One half these ”worst” case values were 
actually used in an attempt to average the Q matrix over the full 
range of expected maneuvers. 
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Estimation of Q (l,l) 

Assume the aircraft can vary speed by (lOO n.m / hr. max.) 
0.027 n.m. /sec. 

Q = E [W ] where W is a random 
signal caused by target maneuvers. 

Taking 1/2 (0.027) implies 

Q(1,1) = (0.013) (0.013) = 0.0169 
Estimation of Q(2,2) 

Assume the aircraft can accelerate at (ig) max. 

Taking 1/2 (ig) ^ .0027 n.m. x 1 sec. 

sec"^ 

implies that Q(2,2) = (.0027) (.0027) = .00001 

Estimation of Q(3,3) 

Assimie the aircraft can make a (l2° / second) max. turn. 

Taking 1/2 (l2°/sec) = 6°/sec. and using fig. 6 

We have: 



chord length = 2r sin 0/2 

or 0 = 2 sin"l ^ = 0.083° 

30 

Hence. 

Q (3,3) = (. 083 ) (. 083 ) = .00692 



Estimates for the remaining terms of Q were obtained in a like 
manner with Q taking the following* form.. 



0.0169 


0 


0 


0 


0 


0 


0 


0.00001 


0 


0 


0 


0 


0 


0 


0 . 0069 


0 


0 


0 


0 


0 


0 


0.25 


0 


0 


0 


0 


0 


0 


0.0011 


0 


0 


0 


0 


0 


0 


0.0011 
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’■'ifl;. t). CHANGt, fl'I RADAR AZIMUTH DUE TO A TURNING TARGET. 
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Simulations were run to optimize the Q values in order to obtain 
the best filter response over the range of 0°, 3°, 6°, and 12° /sec. 
turns with constant speed targets. 



0.00169 0 0 0 0 0 

0 0.0030 0000 

0 0 0.015 000 

Q exp. = 0 0 0 0.10' 0 0 

‘0 0 0 0 0.0001 0 

0 0000 0.00001 

The Qg matrix values where not increased in magnitude as might 
at first be expected but were set equal to Q^, for the following reasons 



(1) The eight second scan time applies to targets of fairly 
long range compared to target ranges anticipated for sector scan, and 
as such the targets are expected to maneuver as little as possible 
in an attempt to close the ship. 



(2) It is intuitively felt that T(ave.) for aircraft maneuvers 
for this type of target might be expected to be 3-5 seconds vice 8 
seconds . 



(3) The effect of maneuvers acting on a spherical coordinate 
system decreases with increasing range. 
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CHAPTER IV 



CONCLUSIONS 

The radar system simulation presented in this paper is felt over- 
all to be a simplified but realistic model of a sophisticated system 
which could be implemented with current ’’state of the art” hardware. 

The basic main program / subroutine organization divides the 
system into easily recognized functions capable of being extended in 
scope, modified in part, or replaced entirely with a new concept, 
while maintaining the integrity of the remainder of the program. 

Any portion of the model can be reworked to meet special requirements 
of a’ particular problem of interest. 

By the addition of new subroutines which might include phenomena 
associated with clutter, target scintilations , and atmospheric atten- 
uation a more realistic detection probability, given a set of in- 
itial conditions, could be achieved. 

Hence the original objective of building a model for simulation 
of a hypothetical radar system as a tool for analysis of a variety 
of radar system problems has been met. 

As an exercise with practical significance, the simulation was 
used to investigate the relative effectiveness of several types of 
radar tracking filters. 

While no pretension of an extensive analysis is implied, the 
author feels that a reasonable and unequivocal comparison of the 
filters can be made from the material presented. 

Based on the ensemble averages of the squared difference between 
predicted and actual positions the Kalman filter obviously provides 
a somewhat better tracking response for all target tracks tested. 
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The tracking ability of both filters appears to be about equal 
for ^^look alike" targets in close proximity. The ability of the 
tracking routine to resolve these targets depends on using relevant 
gate sizes. The gate size in turn is directly related to the accu- 
racy of the radar measurement of range , and angles , and to the extent 
of maneuver capability of the target. The minimum gate size allowable 
without inducing excessive non-correlation due to observation noise 
would thus be essentially the same for both filters. 

The most significant advantage of Kalman Filter tracking appears 
to be in its ability to track maneuvering targets without dropping 
track. This advantage is a result of its ability to automatically 
increase the gate size if a firm track fails to correlate, thus im- 
proving the probability of correlation on the next scan. 

The requirement for peripheral tracking functions such as track 
firmness and quality would surely be reduced in a system using Kalman 
filtering since the (P) matrix can be used as a measure of this type 
of imformation. This feature affects somewhat the disadvantage of 
increased computation time due to the recursive algebraic equations 
required to calculate the G and P matrices. 

Finally, the dependence of the Kalman filter on the Q matrix 
dictates that an efficient and reliable method for estimating Q for 
a particular target is essential. 
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APPENDIX I 



GRAPHICAL OUTPUT 



Graph 










Page 


Alpha-beta filter 


0°/sec turn rate. 
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Kalman filter 


0^/sec turn rate. 
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Alpha-beta filter 


3°/sec turn rate. 
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Kalman filter 


3°/sec turn rate. 
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Alpha-beta filter 


6°/sec turn rate. 








Uo 


Kalman filter 


6°/sec turn rate. 








i+1 


Alpha-beta filter 


12°/sep turn rate. 








h2 


Kalman filter 


12°/sec turn rate. 








1+3 


Alpha-beta filter 


Crossing targets T 


= 


8 


sec . 


1+1+ 


Kalman filter 


Crossing targets T 


= 


8 


sec . 


1+5 


Alpha-beta filter 


Maneuvering target 


T 


= 


8 sec . 


1+6 


Kalman filter 


Maneuvering target 


T 


= 


8 sec . 


1+7 


Alpha-beta filter 


Crossing targets T 


= 


1 


sec . 


1+8 


Kalman filter 


Crossing targets T 


=: 


1 


sec . 


1+9 
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Interpretation and constraints 

(1) The first eight graphs provide a comparison of the pre- 
diction accuracy of each of the two filters investigated. 

The following rules apply: 

(a) Four different target tracks are employed i.e. 0°, 

3"^, 6^ and 12''/sec. 

(b) Target turning rates are initiated at T = 6 sec. 

(c) Each graph represents the results in range (□), 
azimuth (X), and elevation (^) of a "Monte Carlo" ensemble average 
of 100 runs . 

(d) Sample rate = 1 sec. appropriate to Mode 2 operation. 

(e) Alpha-beta filtering is performed in a cartesian coordinate 
reference frame while Kalman filtering is done in spherical coordinates. 

Filtering of the azimuth coordinate, which happens to be changing 
most rapidly, is about equal for both filters. The Kalman filter 
is obviously superior in range and elevation predictions. 

The following symbol table applies for all succeeding graphs. 

true track 

X predicted position target one. 

^ predicted position target two. 

(2) Graphs nine and, ten represent filtering ability for "look 

•f’' 

alike" targets, i.e. similar range, elevation, and azimuth. 

The ability of both filters, supplemented by a set of correlation 
and association rules, to distinguish between similar targets in close 
proximity leaves something to be desired. A reduced revisit rate i.e. 

(l sec vice 8 sec.) has little or no effect, (see graphs 13 and l4). 
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It seems reasonable therefor to assume that a higher resolution radar 
is the only solution to this, problem. 

(3) Graphs ten and eleven indicate the improved tracking response 
of the Kalman filter for maneuvering targets, aided by a variable gate 
size technique. Investigation of the printed output for this run shows 
that the Kalman filter operating on spherical coordinates consistantly 
predicted the position of the turning target within the limits of the 

V 

gate size. 

In contrast the Alpha-beta filter was unable to follow azimuth 
changes with sufficient accuracy to obtain the necessary correlations 
and as a result dropped track excessively. 
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ALPHA-BETA FILTER TURN RATE 0 DEG/SEC. 

X-TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y-AVERAGE ({PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 



-w- 

cr 

2 > 



-w 















-u- 

K 



-96- 



dd 

-96- 



dC 



'M 



-Bd- 









- 86 - 



— a— 

CSl 

a 

a- 

a 

a 



oci 

oq 

ea- 

□a 

□ a 

— Q a- 

□ a 

IS] 

-^- €D — 



&: 






a □ 
□ 






300 



h00 



300 



000 



37 



Qm mo 006 009 012 015 018 021 

KALMAN FILTER TURN RATE 0 DEG/SEC. 

X - TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y - AVG .(( PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 
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ALPHA- BETA FILTER TURN RATE 3 DEG /SEC. 

X-TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y-AVERAGE ((PREDICTED - ACTUAL )**2) Y-SCALE - 0.20 UNITS/INCH. 
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KALMAN FILTER TURN RATE 3 DEG/SEC. 

X - TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y - AVG. ((PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 
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ALPHA-BETA FILTER TURN RATE 6 DEG/SEC. 

X-TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y-AVERAGE ((PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH, 
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KALMAN FILTER TURN RATE 6 DEG/SEC. 

X - TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y - AVG. ((PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 
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ALPHA-BETA FILTER TURN RATE 12 DEG/SEC. 

X-TIME IN INTEGER SECONDS X-SCALE - 3.00 UNITS/INCH. 

Y-AVERAGE ((PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 
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KALMAN FILTER TURN RATE 12 DEG/SEC, 

X - TIME IN INTEGER SECONDS X-SCALE - 3,00 UNITS/INCH. 

Y - AVG. ((PREDICTED - ACTUAL)**2) Y-SCALE - 0.20 UNITS/INCH. 
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ALPHA-BETA FILTER 
X = Y = N.M. 

X SCALE = Y SCALE =3.00 UNITS/INCH. 
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KALMAN FILTER 
X = Y = N.M. 
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APPENDIX II 



LIST OF VARIABLES 

This appendix contains a listing of the variables used in the 



simulations , 


arranged in alphabetical order. 


Variable 


Definition 


ALPHA 


Position smoothing constant of the OC-B tracker. 


ALT(K) 


Altitude of the target. 


ALTl 


A dummy variable. 


AZBEAML 


The leading edge of the radar beam. 


AZBEAMT 


The trailing edge of the radar beam. 


AZDIFF 


Predicted azimuth - true azimuth. 


AZDOT 


Azimuth rate. 


AZ1+(k) 


+ Vi 

Counter for the number of times the K track 

fails to correlate. 


azg(k) 


th 

Azimuth gate size of the track. 


AZ2(I) 


The ensemble average of ELDIFF **2 


AZRAD(K) 


"fch 

The noisy radar azimuth of the ^ target. 


AZI 


A dummy variable. 


AZT(K) 


th 

The true azimuth of the.K^ target. 


AZY 


A dummy variable. 


BETA 


Smoothing constant for velocity in theo-B tracker 


BOXMAX (I) 


The maximum value allowed for a radar return if 
it is to be associated with a particular range box 


BOXMIN(I) 


The minimum >yalue allowed for a radar return if 
it is to be* associated with a particular range box 


C(I,J) 


Observation / track cor;relation matrix. 


cont(k) 


Azimuth of the leading edge of the radar beam on 
initial detection of the target. 


DELTA 


Normalized velocity vector. 
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Variable 


Definition 


DOTT 


A dummy variable. 


DOTX 

DOTY 

DOTZ 


The predicted velocity components of the K target 

as given by the(X-B tracker. 


ELBEAMB 


Bottom extemity of the radar beam. 


ELBEAMT 


Top extemity of the radar beam. 


ELDIFF 


Predicted elevation - true elevation. 


ELDOT 


Elevation angle rate. 


ELERR 


Radar measurement noise added to the true elevation 
angle . - 


ELG(K) 


" 1 " v» 

Elevation gate size of the K ^ track. 


ELRAD(K) 


• 

Noisy radar value of target elevation.- 


ELT(K) 


True target elevation angle. 


ELTEST 


Parameter for distinguishing between range boxes by 
elevation angle. 


FLAG 


Flag is set by MOVAVE if a target is detected. 


G(I,J) 


Gain matrix for the Kalman filter equations. 


GO 


Used by SCAN 360 to compress full scan. 


H(I,J) 


Observability matrix for Kalman filter equations. 


HD(K) 


Heading of the K^^ aircraft in degrees . 


HPHT 


A dummy variable. 


HT(I,J) 


The transpose of the observability matrix. 


lA 


Range box index. 


lAT(l) 


Used as a counter in conjunction with lAZ. 


lAZ(l) 


Radar return signal counter which is compared to 
a threshold value to determine whether or not a 
detection has occured. 


lAZREF 


Reference azimuth of sector scan. 
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Variable 


Definition 


II 


Elevation beam position indicator. 


INDE(l) 


SEE INDEX (I). 


INDEX (I) 


IJIDEX in conjunction with INDE determine when a track 
failSytwo consecutive times, to correlate thus in- 
itiating drop track prodedure. 


ITEL 


The fixed point form of TEL. 


ITELI 


A d\immy variable. 


ITEMP 


The fixed point form of TEMP. 


ITEMPl 


A dummy variable. 


KC(I) 


A flag used to initialize the Kalman filter for a 
particular observation. 


KER 


Flag to indicate an improper matrix* inverse operation 


LIMB 


Lower limit of the sector scan. 


LIML 


Port limit of the sector scan. 


LIMR. 


Starboard limit of the sector scan. 


LIMT 


Upper limit of the sector scan. 


NB 


The number of runs. 


NBT 


Number of discrete times. 


NO 


Counter for the number. of target bins currently 
active. 


NR(I) 


Flag to indicate known targets in sector scan. 


NT 


Number of targets . 


NUMB 


Sector scan counter. 


NUMl 


Dummy variables. 



NUM2 



NUM3 


• 


NUMSCAN • 


The number of 360 degree scans of the radar. 


NUMTRK 


Number of current active tracks. 
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Variable 


Definition 


OISE 


A dummy variable for the noisy radar observations 
(Range, Azimuth or Elevation.) 


p(i,j) 


The covariance of error matrix. 


PAZRAD 


The predicted radar azimuth of the target. 


PAZY 


A dummy variable. 


PD(I) 


A discrete value of probability from the proba- 
bility of detection table. 


PDET 


Probability of dection. 


PELRAD 


Predicted value of the target elevation. 


PHI(I,J) 


The state transition matrix used in the Kalman 
filter. 


PHI1(I,J) 


The state transition matrix for the one second case. 


PHI8(I,J) 


The state transition matrix for the eight second case 


PHIT(I,J) 


Transpose of the PHI matrix. 


PITCH(K) 


Value of the nose attitude of the target in 

degrees . 


PITCHl 


A dummy variable. 


PKl 9(1, J) 


Storage for P matrix associated vith each track. 


PRADRNG 


Predicted radar range of a target. 


PRED 


A dummy variable used in the Alpha-Beta tracker 
for the predicted value of range, azimuth or ele- 
vation. 


PXRAD 

PYRAD 

PZRAD 


Predicted X,Y,Z, components of target position. 


Q(I,J) 


Matrix representing target maneuver capability. 


Q1(I,J) 

Q8(I,J) 


Q may be set equal to QI or Q8 as required by the 
appropriate scan time. 


R(I,J) 


Covariance matrix of the observation noise. 
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Variable 


Definition 


RADDEG 


Constant used to convert radians to degrees. 


RADIFF 


Predicted range - true range. 


RANDOM 


Uniformly distributed random number. 


RANERR 


Error in measuring target range. 


RANGE(I) 


Discrete value of range used to enter the pro- 
bability of detection table. 


radrng(k) 


Noisy target range as given by the radar. 


RADRNGl 


A dummy variable. 


RDOT 


Range rate. 


rg(k) 


"til 

Range gate size for the K track of the cor- 
relation routine. 


S(I) 


A dummy variable used for temporary storage. 


SMOOTH 


The filtered radar observation in range, azimuth 
or elevation. 


SUM 


Temporary value of the radar beam leading edge. 


T 


Time interval for updating targets equal to scan 
time. 


T1 

T2 


Dummy variables . 


TARGRNG 


True value of target range as given by the target 
generator. 


TEL 

TEST 

TESTT 

TNT 


Temporary sector scan beam reference. 


THETA 


Value of target heading in radians with refer- 
ence to the positive X-axis. 


VEL(K) 


th 

Speed of the K target in n.m./hr. 


VMAG2(l) 


Ensemble average of radiff. 



Variable 


Definition 


VNB 


Fixed point form of NB. 


VPHI2(I) 


Ensemble average of azdiff. 


XP(I,J) 


Predicted state vector. 


XPI 


A diiminy variable. 


XRAD(K) 


X coordinate of the target as given by the 

radar . 


XRADl(K) 

XRAD2(K) 


Stored radar track quantity for graphical output. 


XS 


Smooth state vector. 


XT(K) 


th 

The true X coordinate of the K target 


XTl(K) 

XT2(K) 


Stored target track quantity for graphical output 


YRAD(K) 


Radar Y coordinate of the target. 


YRADl(K) 

YRAD2(K) 


Stored radar track quantity for graphical output . 


YT(K) 


The true Y coordinate of the target. 


YTl(K) 

YT2(K) 


Stored target track quantity for graphical output 


Z(I,J) ' 


Measurable quantities input to the Kalman filter 
equations . 


ZRAD(K) 


The radar Z coordinate of the target. 


zt(k) 


The true Z coordinate of the target. 
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Note: Variables associated with the following input/output and 

matrix algebra subroutines, have not been defined., since 
they are essentially dummy variables . 



1. ADD 

2 . PRINTT 

3 . PROD 
h. READD 

5 . RECIP 

6. TRANS 
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APPENDIX III 



LOGIC FLOW DIAGRAMS 



This Appendix contains a series of Logic Flow diagrams as listed 
helow. 



MAIN PROGRAM 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 



TARGEN 

SCAN 360 

DETECT 

MOVAVE 

CORRASS 

ASSOC 

ABFILT 

KALFILT 

SETUP 

TWS 

MONTE 
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FLOW CHArT symbols 



r 





A connector or terminal. 



An offpage connector. 



A predefined process or module /subroutine. 

A more detailed flow chart of this subroutine 
is also included. 




Input /output other than display. 




Decision. 



Processing , annotation . 
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MONTE 
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APPENDIX IV 



PROGRAM LISTING 
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BO 55 J= 1»10 
55 C(IC»J)=0.0 
160 DO 19 J=1,NT 

PRINT 61»J»NR(J) 

61 FORMAT! /»10X»2HJ=. 12 »6HNR(J) = , 12 
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TEST RUN FOR MONTE CARLO AVERAGES 

THIS PROGRAM TESTS THE iCALMAN FILTER FOR AN OBSERVATION TIME OF 
VMAG2(L)= ENSEMBLE AVERAGE OF RADIFF 
VPHI2(L) = ENSEMBLE AVERAGE OF AZDIFF 
AZ2 (L) = ENSEMBLE AVERAGE OF ELDIFF 



B = THE NUMBER OF RUNS 
NBT = NUMBER OF DISCRETE TIMES 
NB=100 
NBT=20 
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